SACLAY-T00-182, 
TSL/ISV-00-0239 
December 2000 



Factorial correlators: angular scaling within 

QCD jets 



R. Peschanski * ft B. Ziaja ** f Q 

CEA, Service de Physique Theorique, CE-Saclay, 
F-91191 Gif-sur-Yvette Cedex, France 

** Department of Theoretical Physics, 
H. Niewodniczanski Institute of Nuclear Physics, 
31-342 Cracow, Poland 

t High Energy Physics, 
Uppsala University, 
P.O. Box 535, S-75121 Uppsala, Sweden 



Abstract: Factorial correlators measure the amount of dynamical correlation in 
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derivation of factorial correlators for a QCD jet described at the double logarithmic 
(DL) accuracy. We obtain a new angular scaling property for properly normalized 
correlators between two solid-angle cells or two rings around the jet axis. Normalized 
QCD factorial correlators scale with the angular distance and are independent of the 
window size. Scaling violations are expected beyond DL approximation, in particular 
from the subjet structure. Experimental tests are feasible, and thus welcome. 
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1 Introduction 



Large multiplicity fluctuations observed in high energy collisions have been already 
studied for many years Q. Advanced methods of data analysis like e. g. factorial 
moment approach [Q, [3], |4[] have been introduced and implemented for the analysis of 
multiplicity patterns. Finally, they led to discovery of intermittency in multiparticle 
production which refers to the scaling of factorial moments with the size of a single bin 
within the analized pattern || |J. 

Many different models have been proposed for the explanation of the effect ||. 
Some of them suggested that an underlying final state multiparticle cascade may be 
responsible for the scaling of particle moments p, ||. Straightforward calculations 
performed for multiplicative random cascading models 0, |3], [7] led to the qualitative 
predictions for the scaling behaviour of factorial moments which were backed afterwards 
by analyses proposed in the framework of the standard theory of strong interactions 
(QCD) ||. Monte-Carlo simulations based on conventional QCD parton cascading 
tend to describe quite well the effect |§ , confirming the relevance of scaling properties 
in QCD parton cascading. However, discrepancies in the precise comparison with 
analytical predictions remain |J. 

So far, both the phenomenological and theoretical investigations of multiplicity pat- 
terns have concentrated mostly on different kinds of particle moments estimated for a 
single bin ||. There are, however, still intriguing questions remaining about the proper- 
ties of correlations between different bins. The observables related to these correlations 
are expected to reflect the presence of large dynamical fluctuations underlying the pat- 
tern stronger than the averaged observables estimated for a single bin. To investigate 
these bin-bin correlations the factorial correlators [[| [| have been introduced. 

Factorial correlators seem to contain some extra information on multiplicity fluc- 
tuations which may be used to complete and adjust the information obtained on that 
from the standard factorial moment analysis ||10|| . Moreover, the present status of ex- 



perimental investigations H, [11] allows one to expect that the comparison of model 



predictions with the real data will be possible soon. It could help to investigate more 
systematically the validity of QCD Monte-Carlo for the description of fluctuations in 
jets and discuss, using a wider set of data, the relevance and problems of analytical 
QCD calculations. 

However, factorial correlators have been studied only in the framework of phe- 
nomenological models |J. The rigorous analysis in the framework of QCD (even at 
leading log orders) has not been performed so far. 
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Figure 1: Example of phase-space cell for the QCD parton cascade. The two-dimensional cone of 
half-opening angle 9 is placed at both solid angle 9q and azimuthal angle 4> taken with respect to the 
main jet axis. 



This paper aims to fill the gap by presenting the analytical derivation of factorial 
correlators performed for the QCD parton cascade || at the double logarithmic (DL) 
accuracy fl2"l . For simplicity we consider only the fixed as case, expecting that it 



gives good qualitative estimation of scaling exponents as it was realized previously for 
factorial moments case ||. The obtained scaling dependence of the correlators on the 
relative distance between the two solid-angle cells recovers the similar result obtained in 
the framework of random cascading a— model 0. This seems to be a kind of universal 
relation. 

The paper is organized as follows. In the next section we introduce factorial cor- 
relators defined for small solid angle cells around the subjet direction. In section 3, 



using the DL generating functional [12], we derive the inclusive two gluon distribution 
which is neccessary to evaluate the correlators. In section 4 the leading contribution 
to factorial correlator is estimated, and found to obey a scaling law similar as in ref. 
for random multiplicative cascading models. In section 5 we discuss briefly the 
modifications which could come from relaxing some of our approximations: including 
running as or energy momentum conservation. Finally, in section 6 we sum up our 
results and present our conclusions, including suggestions for experimental evaluation 
of the normalized factorial correlators and a discussion of the new QCD scaling law 
found in the DL approximation of QCD. 
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2 Factorial correlators in QCD jets 



Normalized factorial moments F q 0, |3|, f|, ||] designed to study multiplicity fluctuations 
in a given phase-space cell of size 8 are defined as : 



(n(n-l)...(n-q + l)) s 
(n) q 5 



where n is particle multiplicity in the phase-space cell, and the average () is made 
over events. Among other types of fluctuations studied using factorial moments, the 
intermittency regime corresponds to moments which scale with the size of phase-space 
cell as : 

F *W~[j) > (2) 
where A denotes the size of whole available phase-space, and <ft q is a positive scaling 
exponent known also as intermittency exponent. 

In order to study correlations between different phase-space cells one introduces 
factorial correlators F qiiQ2 (known also as multivariate factorial moments) 0, || which 
may be regarded as multidimensional extension of moments (]l|). They take the form : 

, , (njn — 1) . . . (ra — gi + 1) \ Sl n(n - 1) . . . (n - q 2 + 1) |<5 2 ) , , 
F qi , q2 {6i,6 2 ) ( n) ^ ( n )g ' ^ 

where Si and S 2 denote the sizes of two separate phase-space cells. Assuming a multi- 
plicative random cascade underlying the particle production, one predicts a recursive 
relation between the scaling exponents for factorial moments and factorial correlators : 



Afl>92 S y <3i+<?2 rqi r<?2> 



(4) 



where <p qiy q 2 is the intermittency exponent defined for doubly normalized factorial corre- 
lators which correspond to factorial correlator ([$]) divided by factorial moments derived 
for Si and S 2 cells respectively : 



F qi , g2 (Si,S 2 ) /A, 
F qi (Si) F q2 (S 2 ) ~ U 12 



"81, 92 

) 



(5) 



where Si 2 is the relative distance between the two phase-space cells Si and S 2 . Note 
the interesting feature that the dependence on the individual phase-space cells Si and 
S 2 disappears from the factorial correlator when they are normalized as in (|5|). 

For the QCD parton cascade M the phase-space cell is more conveniently choosen 
to correspond to a window in emission solid angle 9 (see e.g. Fig. 1). The window size 
6 is then compared to large scale 9q which denotes the jet emission angle. The window 
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may be either one- dimensional ring of aperture 9 placed at the angle 9$ with respect to 
the sphericity axis, or it may be a two-dimensional cone of half-opening angle 9 placed 
at both the solid angle 9q and azimuthal angle taken with respect to the main jet 
axis. 

It was found out that for QCD angular factorial moments || there is a scaling 
relation : 

,QCD 

F q {9) ~ I J) , (6) 

where the intermittency exponent <j)® CD calculated for simplicity in double logarithmic 
approximation with fixed coupling constant as reads : 

0? CD = --7og +D{q-1). (7) 
q 

Number D denotes the window dimension, and it equals D — 1 for the ring, and D = 2 
for the cone. The coefficient 70 is the QCD anomalous dimension for gluon cascade 
which for fixed as equals 7o = 4Nc^ [jl2|| , where Nq denotes the number of colours. 

In order to extend the intermittent analysis by investigating also the possible cor- 
relations between particle flows measured at two different rings around the subjet axis 
(E, #0) let us introduce the angular factorial correlators for the QCD parton cascade, 
defined as follows. Similarly as for the angular factorial moments one would identify 
the large scale (size of whole available phase-space A) with the respective subjet emis- 
sion angle 9q, an d the small scales S\, 62 with the window apertures 61, 9 2 (cf. Fig. 2). 
We will consider parton flows emitted into two rings placed at the separation angle 
#12 with respect to the subjet axis. The ring openings are 9\ and 9 2 respectively (see 
Fig. 2). In using DL approximation framework, we have to assume that the angles are 
small with respect to the subjet direction. More precisely, we will assume that they 
obey the inequalities : 

01, 02 < 012 « 00- (8) 

Henceforth #01 ~ $02, and the relative bin distance 5 12 i n the one-dimensional approxi- 
mation then corresponds to the angular distance 6*12 = #02 — 0oi between the two rings. 
We will discuss the relevance of this DL approximation in our discussion in section 5. 

Having defined angular factorial correlators, we may now estimate them with a 
good accuracy from the convolution of the inclusive two-particle density 
D( 2 \P; E, #0; ki, k 2 , 0i2, 0i, 02) with the respective multiplicity moments in the phase- 
space cells 9i and #2- Using (as for the case of factorial moments ||) their expression in 
the so-called KNO limit proportional up to constants to gith and q 2 th power of mean 
multiplicities N(ki9i), iV(A; 2 02), we obtain : 

F qiq2 (6 ; 0i2,0i, 2 ) oc / — ^ / j^-D^iP; E,9 -, k 1 ,k2,6 1 2,6 1 ,82)N qi (k 1 8 1 )N q2 (k 2 6 2 ) , (9) 
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Figure 2: Phase-space cells for angular factorial correlators. Parton flows are emitted into two rings 
placed at the separation angle 9 12 with respect to the subjet axis. The ring openings are Q\ and 62 
respectively. 



E 




a) 



b) 



Figure 3: Diagrammatic representation of the single particle inclusive density D^(E/kx,6o/9i) 
(Fig. 3a) and the modified inclusive two particle density (E/ki, E/kz, O0/O1, do/02] #12) (Fig. 3b). 

where E denotes the energy of the subjet. 

The mean multiplicity for QCD parton cascade is dependent on an infrared cut-off 
fi, and it reads : 



However, similarly as for the factorial moment case we expect that the cut-off depen- 
dence will disappear after normalization, i. e. when coming to normalized factorial 
correlators (H). The inclusive two-particle density D^{P\ E,9 ; ki, &2, #12, #2) re- 
mains thus the only unknown quantity, necessary to evaluate Its explicit form will 
be derived in the next section. 



N{k9) 



~ e 




(10) 
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9, 




Figure 4: Diagrammatic representation of the convolution Dp' ex \E, Oq) x 

D l - 2 ^(E/ki,E/k2,0Q/8i 1 8o/82;9i2)- It corresponds to the emission of particles from one subjet 
(E,6q) originating from the main jet (P,9p) into two rings of apertures 61,82 with separation angle 
9\2, placed around the subjet axis. 

3 Inclusive two-particle density 



In order to obtain the explicit form of two particle density D (2) (P; E,9 ; fci, fe, #12, #1, #2) 
to insert it into (|9[) we start with a derivation of a related quantity D^ik^k-i) from 
the QCD parton cascade formalism. 

The inclusive two-particle density D^(k 1 ,k 2 ) i.e. the inclusive density to produce 
two particles of momenta k%, k 2 from a parent particle of momentum p may be derived 



in a convenient way from the generating functional for the QCD parton cascade [12]. 
This functional in DL approximation with fixed as takes the form : 

Z p [u] = u(p)eI d3kM ^^^ (11) 

with the initial condition Z p [u] \{ u =iy — 1- The function u(p) is a probing function while 
the factor M p (k) describes the DL probability of emitting a particle of momentum k 
from a primary particle of momentum p. It reads : 

d 3 kM p (k) =^ff^-e p (k), (12) 

where cut-off ^-function 8 p (k) contains phase-space limitations resulting from possible 
angular and energy ordering between parent particle p and child particle k : 

6 p (k) = {p> k, 9 pk < P , kO pk > 11} , (13) 

where \x is the infrared cut-off. 

The inclusive two-particle density is then defined as a functional derivative of Z p [u] : 

= s Jlf u{k2) \ ( ^, (M) 
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which results in a recursive equation for D^(ki, k 2 ) : 

Df> (k 1 ,k 2 ) = D^(k 1 )D^(k 2 ) - SQnb/hVSQnb/kt)) + / d 3 k M p {k) D^{h,k 2 



*2y, 

(15) 

where D^\ki) D^\k 2 ) denote the single particle inclusive densities for particles k\ and 
k 2 respectively, and the integration limits are defined by (|13|) (max(9i, 6 2 ) < 9 < 9 P , 
max(ki, k 2 ) < k < p). 



Taking into account the scaling properties of the DL phase-space measure fll2|) , one 
obtains an equation for the modified inclusive two-particle density 
D( 2 \E /ki, E /k 2 , 9q/9i, 9o/9 2 ;9i 2 ) with a scaling dependence on the relevant variables 
(choosen here to be p = E, the energy of the subjet, and k\ > k 2 , 9i < 9 2 ) : 

D p 2 Hk u k 2 ) ^ D^(E/k 1 ,E/k 2 ,9 Q /9 1 ,9 /9 2 ;9 12 ) = D^(E/k 1 ,9 /9 1 ) D^(E/k 2 ,9 /9 2 ) 

-S(\n(E/k 1 ))S(\n(E/k 2 ))+ 7o 2 / ^ IT D ^ ( k / k ^ k l k ^ 6 I 9 ^ 6 I 9 * M, ( 16 ) 

Jk x k J 9l2 9 

where we have assumed (cf. (|l|) ) that the intermediate emissions represented by the 
homogeneous term of ([15]) do not generate additional particle flow into either 9\ or 
9 2 window (9 12 < 9 < 9 ). The densities D^(E/k u 9 /9 1 ) and D (1 \E/k 2 ,9 /9 2 ) (cf. 
Fig. 3a) denote here single particle inclusive densities for particles k\ and k 2 respec- 
tively. All notations are as in Fig. 2. The density D^(E/k h E/k 2 , 9 /9 1 , 9 /9 2 ; 9 l2 ) (cf. 
Fig. 3b) precisely represents (still in the one-dimensional approximation) the inclusive 
two particle density to obtain two particles of energies ki, k 2 from the subjet (E,9q) 
separated by the relative angle 6*12. The ring apertures are #1, 9 2 respectively. 

The relation between D^(E/ki, E/k 2 , 9 /9i, 9o/9 2 ; 9\ 2 ) and the two particle density 
(|) is the following : 

D^ 2 \P- E, 9 ; fc ls k 2 , 9 12 , 9,, 9 2 ) = D^ ex) (E, 9 ) ■ D^(E/k u E/k 2 , 9 /9 u 9 /9 2 ; 9 12 ),(17) 

where Dp' ex \E, 9 ) is the exclusive single particle density to produce a subjet of energy 
E placed at the opening angle #0 with respect to the main jet P. 

To sum up, convolution (fl7|) represents the emission of particles from one subjet 
(E,9q) originating from the main jet (P,9p) into two rings of apertures 9i,9 2 with 
separation angle 6*12, placed around the subjet axis (see Fig. 4). 

Introducing new variables : 

h h 9 Q 9 12 9 12 

xi = — , w 12 = — , y 2 = In — , tx = In — , t 2 = In — , (18) 

hi K 2 (7i2 "l "2 



we rewrite (]I5]) as : 

D^^/x^wuMM^) = D m {l/x l ,y 2 + t l )D {l \w l2 /x l) y 2 + t 2 
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- <5(lnl/xi)<5(ln«; 12 /xi) (19) 
+ To / — T / dy' 2 DW(l/x r 1 ,wi2,t 1 ,t 2 ,y r 2 ). 

J x\ X± JO 

In order to solve Eq. (|19"1), we transform it into moment space by means of Mellin 
transform : 

D {2 \oj,w 12 ,t u t 2 ,y 2 ) = f 1 dxxxf 1 D i2 \l/ Xl ,w 12 ,t u t 2 ,y 2 ). (20) 



o 



Then we differentiate both sides of ( 19|) with respect to y 2 . Finally, we obtain : 



d 2 
—D {2 \uj,wi 2 ,t 1} t 2 ,y 2 ) = r(u,wi 2 ,ti,t 2 ,y 2 ) + — D {2 \u, w l2 , t l} t 2 , y 2 ), (21) 
dy 2 uj 

where the function r(ui, Wi 2 , t 12 , y 2 ) reads : 

r{u,w l2 ,t u t 2 ,y 2 ) = / dx 1 xr 1 D (1 \l/x 1 ,y 2 + t 1 )D (1 \w 12 /x 1 ,y 2 + t 2 ). (22) 



Equation ([H]) is an ordinary inhomogeneous linear differential equation. Taking into 
account initial conditions defined by (|19|) , its explicit solution takes the form : 



2 Jl 



D (2) (^,u7i2,ti,t 2 ,Z/2) = r(u,wi 2 ,ti,t 2 ,y 2 ) +^-e° m R(u,wi 2 ,t l ,t 2 ,y 2 



00 

2 2 2 

-8{w 12 - 1) e^ 2 - ^ e^ 2 R(u, w 12 , t u t 2 , 0), (23) 

0J 



where function r(oj,Wi 2} ti,t 2} y 2 ) was defined in ([22]), and function R(oj,Wi 2 ,ti } t 2 ,y 2 ) 
denotes the following indefinite integral of r(u, W\ 2 , t u , y 2 ) : 

R(uj, Wi 2 , t u t 2 , y 2 ) = / dyr(u,w 12 ,t 1 ,t 2 ,y)e » m . (24) 



4 The QCD factorial correlators: Derivation 

As a function of the two-particle inclusive density (^), the convolution of and 
multiplicity correlators in ([)]) can be now expressed in explicit form. After a new change 
of variables : 

E h\ 

/i = log— , s 12 = log— (25) 
k\ k 2 

the equation (||) rewritten for the normalized angular correlators 

p (ft I ft ft I ft \ - Fqi 92 (^0^12 ^1^2) , 9fn 

F " W, " W " BpiWHi)^ (26) 
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takes the form : 

_ roo roo r+ioo shj 

dh ds 12 — e" h D&\uj,w 12 ,t x ,t 2 ,y 2 ) 
Jo Jo J-ioo 2m 

X g-9l70(y2+tl+;i) e -9270(y2+*2+;i+S12) ^27) 

where we substituted D^(l/xi,Wi 2 ,ti,t 2 ,y 2 ) by its inverse Mellin representation : 

D®(l/x 1 ,w 13 ,t 1 ,h,V2)= P-.Q/xx)" D®(uj,w 12 ,t u t 2 ,y 2 ). (28) 

J-ioo 2m 



Now let us calculate term by term the contributions to convolution integral fl2~7|) 
coming from the various components of the modified inclusive two particle distribution 
(f23|) denoted (I), (II), (III) and (IV) as follows: 

D^\u,w m t!,t 2 ,y2) = r(u;, Mia, ti,* 2 , Ife) (!) 

2 2 

+ ll e ^R(u,w 12 ,t 1 ,t 2 ,y 2 ) (II) 

UJ 

- 5(w 12 - 1) e~ y2 (III) 

2 2 

_ ^e^i^ w t t 0) (IV), (29) 

UJ 

These terms have the following Mellin representation : 

/+ioo (Jjjj-^ r+ioo ^jjj^ 
77— / 77— ( 30 ) 

-ica 2m J-ioo 2m 

1 



UJ — UJ\ — UJ 2 



exp ( 7o2/2(— + — ) + ^ 2 si 2 + ?2-ti + —h 
V 0J1 uj 2 U)\ uj 2 



(TT ,il ± mx)( , , , [+*>*> fa! [+^duj 2 

(II) — e» m R(uj,w 12 ,ti,t 2 ,y 2 ) = / — — (31) 

UJ J-ioo 2m J-ioo 2m 

1 1 U)\U1 2 



uj-uj 1 -uj 2 UJ - UJ 1 + UJ 2 



UJl UJ 2 UJ\ UJ 2 



x exp ( 7o2/2(— + — ) + ^12 + —tx + ^t 2 
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2 2 
7 0„„ ,./ s 7 0, 



(///) 5(w 12 - 1) e~ V2 = 5(s 12 )e- y \ (32) 



„2 ^2 



Tn /-moo W/j /-+100 flu* 

(IV) ^ e^> flu;, w 12 , h,t 2 ,0)= / / ^ (33) 

UO J-ioo 27U J-ioo 2m 

1 1 ^Cc^ 



- U! - UJ 2 U- Wi + ^ 2 

/ 2 2 2^ 

x exp ( + cj 2 si 2 + ^ti + ^t 2 

\ 00 UJ 1 UJ 2 



The contributions of fl30D , (|3l"l) , (0), ( p3|) to ( ^7|) may be evaluated using the 



multidimensional saddle point approximation. For first term ( pOf ) its convolution (|2"T|) 
takes the form : 

_ /"oo poo p-\-ioo fyjj 

F^v/hJo/O*) ~ dh ds l2 / —e^r^.w^h^) 

y y JO jO J-ioo 27TZ 

X g-9i7o(y2+ti+'i) e -<?27o(y2+*2+'i+si2) _ 



/*oo /*oo /* 

/ dh / C?Si2 / 
Jo JO J- 



+ io ° du f +ioc duoi r+ io ° duj 2 



2m' J-ioa 2m' J-ioo 2m' 



uj — — uj 2 

x exp {-qij {y 2 +h + h) - q 2 'j (y 2 + t 2 + h + s 12 ) + u h) 

X exp ( 7^2 (— + — ) + ^2512 + —ti + ^t 2 ) . 



2 721 . . 

U\ LU 2 LOi UJ 2 



(34) 



After performing the integral over the u;-pole in ( |3"o] ) one obtains : 



roc poc p+ioa ^ p+ioo ^ 

F L q?(0o/0i, 60/62) ~ / dh \ ds 12 / — — / —^exp(S(uJi,U2,h,si2,ti,t2,y2)) (35) 



which in the saddle point approximation may be estimated as : 

r-oo r-oo r+ioo d(jji r+ioo duJ 2 

/ dl x / ds 12 / — - / -—exp(S{uj 1 ,u 2 ,l 1 ,s 12 ,t 1 ,t 2 ,y 2 )) 
Jo Jo J-ioo 2m J-ioo 2m 

exp ( S(u>i, o; 2 , ii, Si 2 ,ti,t2) 2/2) |^.=o -^-=0 -&§-=o -&§_=o ) 

\ dli ' 3s 12 ' du>\ ' 8^2 / 
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x det 



d 2 S 


d 2 S 


d 2 S 


a 2 s 




8/f 


dhds\2 


dhduji 


dl\du>2 


d 2 S 


d 2 S 


d 2 S 


d 2 S 










ds 12 du>2 




d 2 S 


a 2 s 


d 2 S 


d 2 S 




dujidli 




dtj'f 


du>idu)2 




d 2 S 


d 2 S 


d 2 S 


d 2 S 


) 


duJ2dli 


du)2&S\2 







-1/2 



(36) 



|S =0 , s^ =0 , |« =0, #^=0 



For the function S t (uji, u 2 , 1%, Si 2 , t±, t 2 , y 2 ) reads : 
S J (a>i,a>2,Zi,si2,*i, £2,2/2) = -qilo^ + h + k) - q2lo(y2 + h + h + s 12 ) (37) 



1 1 



T6\ 



+ (wi + cu 2 )/i + 7oya(— + — ) + ^ 2 si2 + — ti + — 1 2 . 



LO-l L0 2 



UJ 2 



Hence after evaluating ( 35|) one obtains : 

1 



7o 



To 



To 



To 



4tt 2 



exp { 1/2 To 9i H To 92 +*i To 9i + *2 To 92 ? • (3 



,9i 



92 



.91 



.92 



For the second term (|31|), convolution fl2"T|) takes the form : 



^ 2 (0oMAM) ~ dh ds 12 , 

y y JO JO J-ioo Z7U 



x g-gi7o(y2+*l+'l) e -9270(y2+*2+'l+Sl2) 



+lOO W/.i ^2 7 2 



too r+ioo (Jjjj r+ioo ^ f+ io ° duj 2 

dli / dsi 2 / ; / r / : 

jo J-ioo 27ri J-ioo 2iti J-ioo 2ni 

1 1 0JiUJ 2 



Ct>l + UJ 2 



10 — uji — uj 2 uj 
x exp (-9170(2/2 + *i + Zi) - <?2To(?/2 + *2 + Zi + s 12 ) + w Zi) 

x exp ( To2/2(— + — ) + ^12 + — h + ^ 2 J • 

V OJi UJ 2 LOi UJ 2 I 



(39) 



Since there are two u>-poles in ( |39|) the integration over w will give rise to two separate 
saddle point exponents : 



-f-zoo 



du 



2iri oj — uj\ — L02 td — 



UJ1UJ2 



ujiLu 2 



which have to be evaluated separately. One obtains : 

- ioc dui f +iQC duj 2 



Ftjeo/^W- I dh ds X2 , 

JO JO J —too Z,tt L 



(40) 



(41) 



[exp( S IIa (uJi,u> 2 , Zi, S12, h, t 2 , y 2 )) - exp( S nb (u h u 2 , h, s 12 , t 1} t 2 , y 2 ))} 
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where S IIa , S IIb read : 

S IIa {u u u 2 , h, s 12 , h, t 2 , y 2 ) = -<?i7o(2/2 + h + h) ~ <?27o(?/2 + h + h + s 12 ) 

+ (uji + uj 2 )h +Joy 2 { 1 ) +uj 2 s 12 

0J1 uj 2 

+ i tl+ 2i t2+ i n ^ „ ,42) 

tui cj 2 ^1 + ^1^2 + ^2 



5 (^l,^2,/l,Sl2,tl,t2,l/2) = -<?1 7o(?/2 + tl + h) ~ q2lo(V2 + h + h + S 12 ) 

^1^2 , 2/1,1 



/1 + 7o2/2( 1 ) +CU2S12 



C^l + ^2 k>l W 2 

70, TO. , I ^1^2 ^ 

+ — ti + — t 2 + ln^ (43) 

Co>l Co>2 ^1 + 0OiUJ 2 + UJ 2 



Hence, after evaluating (|39| ) one obtains : 



F I'Vo/e^BM- (44) 



9l 92 

9i 92 



4tt 2 g? + gig 2 + gl 



exp 



,7o ( 7o \ . , / 7o \ , , / 7o 

2/2 7o gi H 7o 92 + h[ 70 q\ + 1 2 [ 70 92 

qi 92 / \qi J \g2 



(gl +g2)g2 r / 7o , . A , , ( 7o gl \ . , (lo 
2 ex P fa i 7o (gi + qi) ) + h[ -70 gi 7 ; r ) + t 2 [ 70 g2 

qf \gl+g2 / V g2(gl+g2) \g2 



Similarly, for the third fl32|) and fourth term (^), their convolutions fl2*7|) result in : 



7n 2 



q q Jo Jo J-ioo 2m 



-+ioo duo poo r+ioc dp 

dh I 77— I ds 12 



J-ioo 2m' Jo J-ioo 2ni 

x exp (-9170(3/2 + h + h) - q 2 'j (y 2 + t 2 + h + s 12 ) + u h) 

x exp ( ^-y 2 +ps l2 \ , (45) 



_ poo poo p+ioo dcj 1 fn 

(e /e u e /6 2 ) ~ dh ds l2 / e <*i° e -^ R{ui,w 12 ,t u t 2 ,Q) 

' JO JO J— ioo ZlXl LO 



x e -HlQ{y2+tl+h) e ~927o(y2+t2+'l+Sl2) 



00 poo 



r+ioo dio r+i°° dux r+ io ° du! 2 
dh / ds 12 / -— 



JO J-ioo 2ni J-ioo 27Ci J-ioo 2m 
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1 1 C^iCg>2 

u - u x - U) 2 U- Wi + U) 2 

x exp (-gi7o(j/ 2 + *i + ii) - <?27o(l/2 + *2 + h + s i2 ) + w Zi) 

( 2 2 2 \ 

uo Ui uj 2 J 

(46) 

The respective saddle functions obtained similarly as in fl4~0|) and (f4i~D are the following : 

S IH (u,p, li,s J 2,t 1 ,t2,y 2 ) = -qilo(y2 + ti + k) - q2lo(y2 + t 2 + k + s 12 ) 

+ u)l x + —V2+psi 2 , (47) 

to 1 

S Wa (u 1 , U 2 , l h S12, tl, *2, 1/2) = -gi7o(?/2 + ti + /1) - g 2 7o(l/2 + ^2 + ^1 + S12) 

+ (lUi + u 2 )h + 7o?/2 ; hW 2 Si2 

+ C^2 

+ i fl+ 2i t2+ i n ">■» ,48) 

Wi c<j 2 u;f + WxCf 2 + ^ 2 

S ,ZV6 (a'i, u 2 , h, S12, ti, £2, 1/2) = -gi7o(?/2 + *i + ii) - <?27o(2/2 + *2 + ^1 + Sl 2 ) 

^1^2 2/1 1 

+ ; £i + 7o2/2( 1 ) +^2Si 2 

u)\ + uo 2 ojx u 2 

+ — *i H t 2 + ln^- ■ — 2- ( 49 ) 

Wl CJ 2 Ujf + CJi^ + ^2 



Hence after evaluating fl4"5|), (|46D one obtains : 



F^ Q2 (0o/6i,0o/02) ^ (51) 

exp 



91 92 ' 

_1 9192 

4?r 2 <? 2 + gig 2 + gf 
(91 + 92)92 



V2 [ — ^7 7o (91 + 92) J + *i ( — - 7o 9i ) + h ( — - 7o 92 

91+92 / V 9i / V 92 



— , exp 

9? 



2/2 ( — - 7o (91 + 92) ) + h ( -70 91 7-^-7 — r ) +^2 ( — - 7o 92 

91+92 / V 92(91+92)/ \92 , 



We now have to sum up contributions (0), (|44]) , (|50D, ([51]) according to (|27|) , (j2"9|). 
Hence, finally, normalized angular correlators (^6]) read : 

(52) 



iU(MMo/fc) - -LAe y2 (^- 70(9l+92) ) +<1 (^---) +t2 (^---) 

47T Z 



1 1/2 f ^-f^ -7o (91+92)) -t\ (70 9i) -*a (70 92) 
4^ e 

1 (A + 1) e'' 2 (^" 7 ° 91+ ^ 70 92 ) +tl fe~ 70 91 ) + * 2 fe~ 70 92 ) . 



471-2 
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After dividing d52|) , in analogy to @, by the product F qi (9 /9i) F q2 (9o/9 2 ) one 
obtains : 

Wgo/MWgg) J_ , /^\^ + --^-^ _ j_ 

F qi {9 Q /9 x )F q2 {9 Q /9 2 ) ~ 4tt 2 l^J 4^ + J 

+ *?UJ UJ • (53) 

where 9l+(?2 = 70/ (<?i + g 2 ) - 7o(<7i + 92) and gi = 70/gi - 70 gi, (j) q2 = 7o/<?2 - 7o <?2- 
Note that the quantity <f> qi + q2 — 9l — 92 has been denoted <fi qi;q2 in formula (Q). 



5 Angular scaling of factorial correlators 

In order to discuss the physical properties of QCD factorial correlators, let us rewrite 
formula (|52|) in a different form 




• (54) 



4vr 2 

The physical interpretation of the two brackets contributing to formula ([54] ) is quite 



simple. Considering the first one which is dominant at large values of ^ , ^ (§) , it 
corresponds to the contribution coming from the full development of the parton cascade 
(minus the value when 9\ 2 = #o, i-e. substracting the effect of cascading before 9q). 
Indeed, due to the QCD constraints of angular ordering, the angular ordered path 
#0 —> #12 is populated by fluctuations with order q t + q 2 , while the remaining separated 
paths from # 12 — > 9i and #12 — > #2 corresponds to the individual fluctuation patterns 
with order q\ and q 2 . In QCD framework at DLA, this contribution is similar to the 
behaviour of the random cascading models. 



It is clear that this first term in equation fl54j ) implies specific angular scaling proper- 
ties of QCD jets (at DLA). Normalizing this term by the product F qi (9o/9i) F q2 (9 /9 2 ) 
gives : 



Fgi g 2 (^o/^l; ^0/^2) 

F qi (9 /9 1 )F q2 (9 /9 2 ) 




(55) 
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The scaling properties of flog) can be expressed by the following three items 

i) The normalized correlator (at DLA) depends only on the angular separation # 12 , 
and thus is independent of the window sizes 9i,9 2 , 

ii) It obeys a scaling law as a function of the ratio J^-, 

hi) The scaling exponent is related to the ones of the factorial moments by <fi qim = 

091+92 ~~ 091 ~ 092' 

Such a prediction is similar to the one of random cascading models, which has 
previously (|T3|] been discussed for soft hadronic mult ipro duct ion. In that case, the 
property i) has been verified, while the dependence ii) showed some bending and iii) 
was largely violated since from observation in some range of resolution it appeared 
that 4> qi , q2 3> 4> qi +q 2 — 09i — 092- We know now that the multiplicity fluctuations in 
soft hadronic mult ipro duct ion are influenced by Bose-Einstein enhancements. It would 
thus be interesting to measure by comparison the normalized correlators in jets, where 
the dynamics is more directly related to perturbative (and resummed) QCD properties. 
The experimental analysis can be done as an extension of what was done for angular 
factorial moments 0, where window rings around the jet axis have been considered as 
phase-space slices. 



Interestingly enough, a second contribution appears in formula ( |i>iD which also has 
a simple physical interpretation. On the contrary with the first term, the exponents 
09i —Jo/<li an d <p qi —70/91 mean that the parton cascading structure during the second 
step of the process related to the separated paths from 6*12 — > Q\ and 9i 2 — > 2 is damped, 
since the corresponding fractal dimensions 70/91 and 70/91 are cancelled from the 
intermittency exponents. This corresponds to the probability of having a contribution 
of parton jets directly into the windows of observation. This contribution is obviously 
subdominant at DLA, since the exponents are smaller. It would lead to a violation of 
the scaling properties i)— iii) . 

However, it remains to be found whether, beyond the DLA approximation, such a 
contribution could be in practice larger than the first one. In particular, the lack of 
DLA exponentiation could be compensated by the strong contribution to multiplicity 
fluctuations of subjets directly hitting the observation windows. This study is beyond 
the scope of our paper devoted to the analysis of the DLA approximation but is deserved 
in the future. 
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6 Summary 



We presented the analytical derivation of factorial correlators performed for the QCD 
parton cascade at the double logarithmic (DL) accuracy. For simplicity we considered 
only the fixed as case, expecting that it gives good qualitative estimation of scaling 
exponents as it was realized previously for factorial moments. The scaling dependence 
of the correlators on the relative distance between the two solid-angle cells recovered 
the similar result obtained in the framework of random cascading a— model 0, [7]], and 
seems to be a kind of universal relation. 

However, it remains to be found whether the scaling holds also beyond the DLA 
approximation, where the contribution to multiplicity fluctuations coming from subjets 
directly hitting the observation windows may be dominant. This study is beyond the 
scope of our paper devoted to the analysis of the DLA approximation but is deserved 
in the future. 

It would also be useful to compare these predictions with QCD Monte-Carlo sim- 
ulations (based on parton showers). It is already known that there is a noticeable 
difference between QCD predictions at DLA and QCD Monte-Carlo simulations for 
factorial moments, these predictions being in better agreement (but not perfect) with 
data. Since the origin of this discrepancy is not well understood at present, the study 
of factorial correlators could be useful for identifying the problem. 
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